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Abstract 

We report substantial progress in the study of separability functions and their application to 
the computation of separability probabilities for the real, complex and quaternionic qubit-qubit 
and qubit-qutrit systems. We expand our recent work (J. Phys. A 39, 913 [2006]), in which the 
Dyson indices of random matrix theory played an essential role, to include the use of not only the 
volume element of the Hilbert-Schmidt (HS) metric, but also that of the Bures (minimal monotone) 
metric as measures over these finite-dimensional quantum systems. Further, we now employ the 
Euler-angle parameterization of density matrices (p), in addition to the Bloore parameterization. 
The Euler-angle separability function for the minimally degenerate complex two-qubit states is 
well-fitted by the sixth-power of the participation ratio, R(p) = ^ 2 . Additionally, replacing 
R(p) by a simple linear transformation of the Verstraete-Audenaert-De Moor function (Phys. Rev. 
A, 64, 012316 [2001]), we find close adherence to Dyson-index behavior for the real and complex 
(nondegenerate) two-qubit scenarios. Several of the analyses reported help to fortify our conjectures 
that the BS and Bures separability probabilities of the complex two-qubit states are ^ « 0, 242424 
and 1680 (y~ 1 ) ~ 0.733389, respectively. Employing certain regularized beta functions in the role of 
Euler-angle separability functions, we closely reproduce-consistently with the Dyson-index ansatz- 
several BS two-qubit separability probability conjectures. 
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I. INTRODUCTION 

For several years now, elaborating upon an idea proposed in [1], we have been pursuing 
the problem of deriving (hypothetically exact) formulas for the proportion of states of qubit- 
qubit [2] and qubit-qutrit [3] systems that are separable (classically-correlated) in nature [4- 
11]. Of course, any such proportions will critically depend upon the measure that is placed 
upon the quantum systems. In particular, we have — in analogy to (classical) Bayesian 
analyses, in which the volume element of the Fisher information metric for a parameterized 
family of probability distributions is utilized as a measure ( "Jeffreys' prior" ) [12] — principally 
employed the volume elements of the well-studied (Euclidean, flat) Hilbert-Schmidt (HS) and 
Bures (minimal monotone or symmetric- logarithmic-derivative [SLD]) metrics (as well as a 
number of other [non-minimal] monotone metrics [8]). 

Zyczkowski and Sommers [13, 14] have, using methods of random matrix theory [15] (in 
particular, the Laguerre ensemble), obtained formulas, general for all n, for the HS and 
Bures total volumes (and hyperareas) of n x n (real and complex) quantum systems. Up 
to normalization factors, the HS total volume formulas were also found by Andai [16], in 
a rather different analytical framework, using a number of (spherical and beta) integral 
identities and positivity (Sylvester) conditions. (He also obtained formulas — general for any 
monotone metric [including the Bures] — for the volume of one-qubit [n = 2] states [16, sec. 
4]-) 

Additionally, Andai did specifically study the HS quaternionic case. He derived the HS 
total volume for n x n quaternionic systems [16, p. 13646], 

v hs = (2n-2)!7r n ~ n n n-2/ 2 -x, (1) 

giving us for the two-qubit (n = 4) case that will be our specific initial interest here, the 
2 7- dimensional volume, 

7T 12 1 7T 12 

" . ~ 2 93352 • 10~ 18 (2) 

7776000 40518448303132800 315071454005160652800000 ~ ' y ! 

(In the analytical setting employed by Zyczkowski and Sommers [13], this volume would 



appear as 2 12 times as large [16, p. 13647].) 
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If one then possessed a companion volume formula for the separable subset, one could 
immediately compute the HS two-qubit quaternionic separability probability (P^t) by taking 
the ratio of the two volumes. In fact, following a convenient paradigm we have developed, 
and will employ several times below, in varying contexts, we will compute P^t as the 
product (R1R2) of two ratios, R\ and R2. The first (24-dimensional) factor on the left-hand 
side of (2) will serve as the denominator of R\ and the second (3-dimensional) factor, as the 
denominator of R 2 - The determinations of the numerators of such pairs of complementary 
ratios will constitute, in essence, our (initial) principal computational challenges. 

A. Bloore parameterization of density matrices 

One analytical approach to the separable volume/probability question that has recently 
proved to be productive [17] — particularly, in the case of the Hilbert-Schmidt (HS) metric 
(cf. [18]) — makes fundamental use of a (quite elementary) form of density matrix parame- 
terization first proposed by Bloore [19]. This methodology can be seen to be strongly related 
to the very common and long-standing use of correlation matrices in statistics and its many 
fields of application [20-22]. (Correlation matrices can be obtained by standardizing covari- 
ance matrices. Density matrices have been viewed as covariance matrices of multivariate 
normal [Gaussian] distributions [23] . Covariance matrices for certain observables have been 
used to study the separability of finite-dimensional quantum systems [24]. The possible 
states of polarization of a two-photon system are describable by six Stokes parameters and 
a 3 x 3 "polarization correlation" matrix [25].) 

In the Bloore (off-diagonal scaling) parameterization, one simply represents an off- 
diagonal ij-entry of a density matrix p, as pij = y/puPjjWij, where Wij might be real, complex 
or quaternionic [26-28] in nature. The particular attraction of the Bloore scheme, in terms 
of the separability problem in which we are interested, is that one can (in the two-qubit 
case) implement the well-known Peres-Horodecki separability (positive-partial-transpose) 
test [29, 30] using only the ratio, 



rather than the four (three independent) diagonal entries of p individually [11, eq. (7)] [17, 




(3) 



eq. (5)]. 
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Utilizing the Bloore parameterization, we have, accordingly, been able to reduce the prob- 
lem of computing the desired HS volumes of two-qubit separable states to the computations 
of one-dimensional integrals (33) over /i e [0, oo]. The associated integrands are the products 
of two functions, one a readily determined jacobian function (corresponding, first, to 

the transformation to the Bloore variables Wij and, then, to n) and the other, the more 
problematical (what we have termed) separability function S HS '(/i) [11, eqs. (8), (9)]. 

In the qnbit-qutrit case [sec. II C], two ratios, 

PllP55 P22P66 /.X 

v\ = , v 2 = , (4) 

P22P44 P33P55 

are required to express the separability conditions (choosing to compute the partial transpose 
by transposing four 3x3 blocks, rather than nine 2x2 blocks of the 6x6 density matrices),, 
but analytically the corresponding HS separability functions also appear to be univariate in 
nature, being simply functions of either v x or of u 2 singly, or the product [17, sec. Ill], 

7] = v x v 2 = • (5) 

P33P44 

B. Euler-angle parameterization of density matrices 

Here, one can again divide the set of parameters into two groups, in a natural manner 
(that is, the diagonal and off-diagonal parameters in the Bloore framework). Now, the 
two sets are composed of the eigenvalues of p and of the Euler angles parameterizing the 
associated unitary matrix of eigenvectors [3, 31]. 

With both forms of parameterizations we have discussed, one can obtain the total volume 
of n x n quantum systems as the product of integrals over the two complementary sets [13, 14] . 
But this direct approach no longer holds in terms of computing the separable volume. So, 
we have evolved the following general strategy [11, 17]. We integrate over the larger set 
(off-diagonal or Euler-angle parameters), while enforcing separability conditions, leaving us 
with separability functions that are functions of only the smaller set of parameters (diagonal 
entries or eigenvalues). Doing so, of course, substantially reduces the dimensionality of the 
problem. 

We are, then, left with such separability functions and the ensuing task of appropriately 
integrating these functions over the remaining parameters (diagonal entries or eigenvalues), 
so as to obtain the requisite separable volumes. 
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C. Immediately preceding studies 



In our extensive numerical (quasi-Monte Carlo integration) investigation [11] of the 
9-dimensional and 15-dimensional convex sets of real and complex 4x4 density matri- 
ces, we had formulated ansatze for the two associated separability functions and 
^compiexif 1 )) > proposing that they were proportional to certain (independent) incomplete beta 
functions [32], 

CP? 

B^(a,b)= ^{l-uf-^du, (6) 
Jo 

for particular values of a and b. However, in the subsequent study [17], we were led to 
somewhat modify these ansatze, in light of multitudinous exact /ewer-dimensional results 
obtained there. Since these further results clearly manifested patterns fully consistent with 
the Dyson index ("repulsion exponent") pattern {j3 = 1,2,4) of random matrix theory 
[33], we proposed that, in the (full 9-dimensional) real case, the separability function was 
proportional to a specific incomplete beta function (a — |, b — 2), 

SZSiM « B^ 1 -, 2) = ^(3 - ^ 2 )/i = 7^(^33 - PiiP 44 ) (7) 

and in the complex case, proportional, not just to an independent function, but simply to the 
square of S^(fi). (These proposals are strongly consistent [17, Fig. 4] with the numerical 
results generated in [11].) This chain of reasoning, then, immediately suggests the further 
proposition that the separability function in the quaternionic case is exactly proportional to 
the fourth power of that for the real case (and, obviously, the square of that for the complex 
case). It is that specific proposition we will, first, seek to evaluate here. 



D. Objectives 

We seek below (sec. II A) to further test the validity of our Dyson-index ansatz, first 
advanced in [17], as well as possibly develop an enlarged perspective on the still not yet fully 
resolved problem of the two-qubit HS separability probabilities in all three (real, complex and 
quaternionic) cases. (In [17], we proposed, combining numerical and theoretical arguments- 
not fully rising to the level of a formal demonstration-that in the real two-qubit case, the 
HS separability probability is ^, and in the complex two-qubit case, J|.) A supplementary 
treatment of the truncated quaternionic scenario {(5 = 3) is presented in sec. II B. 
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We invesigate related separability-function questions in the qahit-qutrit framework, again 
making use of the Hilbert-Schmidt metric (sec. II C), and, also in the two-qubit setting, em- 
ploying the Bures (minimal monotone) metric (sec. II D). Since it becomes more problem- 
atical to obtain separability functions in the Bures case, we explore-as originally proposed 
in [34]-the use of the SU(4) Euler-angle parameterization of Tilma, Byrd and Sudarshan 
[31] for similarly- minded purposes (sec. III). 

II. BLOORE- PARAMETERIZATION SEPARABILITY FUNCTIONS 
A. Quaternionic two-qubit Hilbert-Schmidt analysis 

Due to the "curse of dimensionality" [35, 36], we must anticipate that for the same 
number of sample ("low-discrepancy" Tezuka-Faure (TF) [37, 38]) points generated in the 
quasi-Monte Carlo integration procedure employed in [11] and here, our numerical estimates 
of the quaternionic separability function will be less precise than the estimates were for the 
complex, and a fortiori, real cases. (An interesting, sophisticated alternative approach to 
computing the Euclidean volume of convex bodies involves a variant of simulated annealing 
[39] (cf. [40]), and allows one — unlike the Tezuka-Faure approach, we have so far employed— 
to establish confidence intervals for estimates.) 

Our first extensive numerical analysis here involved the generation of sixty-four million 24- 
dimensional Tezuka-Faure points, all situated in the 24-dimensional unit hypercube [0, l] 24 . 
(The three independent diagonal entries of the density matrix p — being incorporated into 
the jacobian J(n) — are irrelevant at this stage of the calculations of S^ t (fi). The 24 [off- 
diagonal] Bloore variables had been linearly transformed so that each ranged over the unit 
interval [0,1]. The computations were done over several weeks, using compiled Mathematica 
code, on a MacMini workstation.) 

Of the sixty-four million sample points generated, 7,583,161, approximately 12%, corre- 
sponded to possible 4x4 quaternionic density matrices — satisfying nonnegativity require- 
ments. For each of these feasible points, we evaluated whether or not the Peres-Horodecki 
positive-partial-transpose separability test was satisfied for 2,001 equally-spaced values of 

[0,1]. 

Here, we encounter another computational "curse" , in addition to that already mentioned 
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pertaining to the high-dimensionality of our problem, and also the infeasibility of most (88%) 
of the sampled Tezuka-Faure points. In the standard manner [15, eq. (5.1.4)] [27, p. 495] 
[41, eq. (17)] [42, sec. II], making use of the Pauli matrices, we transform the 4x4 
quaternionic density matrices — and their partial transposes — into 8x8 density matrices 
with [only] complex entries. Therefore, given a feasible 24-dimensional point, we have to 
check for each of the 2,001 values of /i, an 8 x 8 matrix for nonnegativity, rather than a 
4x4 one, as was done in both the real and complex two-qubit cases. In all three of these 
cases, we found that it would be incorrect to simply assume — which would, of course, speed 
computations — that if the separability test is passed for a certain /i , it will also be passed 
for all n lying between po and 1. This phenomenon reflects the intricate (quartic both in p 
and in the Bloore variables Wi/s, in the real and complex cases) nature of the polynomial 
separability constraints [11, eq. (7)] [17, eq. (5)]. 



1. Estimated separability function and probability 

In Fig. 1 we show the estimate we, thus, were able to obtain of the two-qubit quaternionic 
separability function <S>£fi(/i), in its normalized form. (Around /i — 1, one must have the 
evident symmetrical relation S HS (fi) = S HS (j i ).) Accompanying our estimate in the plot 
is the (well-fitting) hypothetical true form (according with our Dyson-index ansatz [17]) of 
the HS two-qubit separability function, that is, the fourth power, ^(3 — fJ 2 )/^ , of the 
normalized form of iS^f z (p). 

For the specific, important value of /i = 1-implying that pupu = p22P33~the ratio (i?i) 
of the 24-dimensional HS measure (m sep = R™ umer } assigned in our estimation procedure to 
separable density matrices to the (known) total 24-dimensional HS measure (m to t = Rf enom ) 
allotted to all (separable and nonseparable) density matrices is R± = 0.123328. The exact 
value of m sep is, of course, to begin here, unknown, being a principal desideratum of our 
investigation. On the other hand, we can directly deduce that m tot = Rf enom = 77 ^ Q00 ~ 
0.118862 — our sample estimate being 0.115845 — by dividing the two-qubit HS quaternionic 
2 7- dimensional volume (2) obtained by Andai [16] by 

R " en ° m = 2 I J ^ d » - - 40518448303132800 * 2 - 46801 ' 10 " " * = * 

(8) 



8 




0.4 0.6 



FIG. 1: Estimate — based on 64,000,000 sampled 24-dimensional points — of the normalized form of 
the two-qubit quaternionic separability function SS^ t (fi), along with its (well- fitting) hypothetical 
true form, the fourth power of the normalized form of <S^^(//), that is, f|(3 — p?)p\ 




FIG. 2: The univariate quaternionic jacobian function J qU at{p) 



Here, 3^ qna t{n) is the quaternionic jacobian function (Fig. 2), obtained by transforming the 
quaternionic Bloore jacobian ^11^22^33(1 — P11 — P22 — P33)) 2 , /? = 4, to the /x variable by 
replacing, say P33 by /i, and integrating out pn and p22- (We had presented plots of J re ai(p) 
and Jcom P iex{.p) in [11, Figs. 1, 2], and observed apparently highly oscillatory behavior 
in both functions in the vicinity of fi — 1. However, a referee of [17] informed us that 
this was simply an artifact of using standard machine precision, and that with sufficiently 
enhanced precision [only recently available for plotting purposes in Mathematica 6.0]-as now 
employed in Fig. 2-the oscillations could be seen to be, in fact, illusory.) We theoretically 
can obtain the two-qubit quaternionic separability probability P^ e p/ quat by multiplying the 
true value (which we do not beforehand know, but seek) of the ratio R\ by a second (known, 
computable) ratio R 2 - The denominator of i? 2 has already been given (8). The numerator 
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of R 2 is the specific value 

R™ mer = 2 / J aua t(u)( -(3- u 2 )/i) du = « 1.67128 • 1(T 17 , (9) 

2 u qua t\wy 2 \ vjvj v 358 347086242825680000 ,KJ 

where, to obtain the integrand, we have multiplied (in line with our basic [Bloore- 
parameterization] approach to the separability probability question) the quaternionic ja- 
cobian function (Fig. 2) by the (normalized) putative form of the two-qubit quaternionic 
separability function. (Note the use of the f3 — 4 exponent.) 

The counterpart of R% umer in the 9-dimensional real case is 151 2 00 and in the 15- 
dimensional complex case, ggggp^ooo ■ We further note, regarding the last denominator, 
that 

I ii \ r(i6) 

99891792000 = -^-f (10) 

V 2 J r (") 

is the coefficient of /i 2 in ll!L^ 1 (//) and 151 2 200 = 75600 plays the exact same role in 6!Lg(/i), 
where is a generalized (a = 4) Laguerre polynomial (see sequences A062260 and 

A062140 in the The On-Line Encylopaedia of Integer Sequences). (Also, as regards the de- 
nominator of (9), ^SS 8000 " = 99891792000.) Zyczkowski and Sommers had made 
use of the Laguerre ensemble in deriving the HS and Bures volumes and hyperareas of 
n-level quantum systems [13, 14]. Generalized (associated/Sonine) Laguerre polynomials 
["Laguerre functions"] have been employed in another important quantum-information con- 
text, in proofs of Page's conjecture on the average entropy of a subsystem [43, 44].) 
We, thus, have, for our two-qubit quaternionic case, that 



R numer 125769 
idenom 



R 2 = —, = ~ 0.677179. (11) 

R denom \%^2h 



(The real counterpart of R 2 is ~ 0.76854, and the complex one, || ss 0.717172. Ad- 
ditionally, we computed that the corresponding "truncated" quaternionic [45] ratio — when 
one of the four quaternionic parameters is set to zero, that is the Dyson-index case (3 = 3— 
* s i06376244975^ 2 ~ 0-692379. Thus, we see that these four important ratios R 2 {(3) monotoni- 
cally decrease as (5 increases, and also, significantly, that the two ratios for odd values of (5 
differ qualitatively — both having tt 2 in their denominators — from those two for even (3.) 

Our quasi-Monte Carlo (preliminary) estimate of the two-qubit quaternionic separability 
probability is, then, 

P s H ep/ g uat ~ R1R2 = 0.0813594. (12) 
10 



Multiplying the total volume of the 27-dimensional convex set of two-qubit quaternionic 
states, given in the framework of Andai [16] by (2), by this result (12), we obtain the 
two-qubit quaternionic separable volume estimate V^, quat « 2.38775 • 1CT 19 . 

Our 24-dimensional quasi-Monte Carlo integration procedure leads to a derived estimate 
of (the total 27-dimensional volume) V q ^ t , that was somewhat smaller, 2.85906 • 10~ 18 , than 
the actual value 2.93352- 1CT 18 given by (2). Although rather satisfying, this was sufficiently 
imprecise to discourage us from further attempting to "guestimate" the (all-important) 
constant (Ri) by which to multiply the putative normalized form, (|(3 — /i 2 )/i) 4 , of the 
quaternionic separability function in (9) in order to yield the true separable volume. In our 
previous study [17, sec. IX. A], we presented certain plausibility arguments to the effect that 
the corresponding R\ constant in the 9-dimensional real case was 1 | 1 5 ^ 2 = (^f-)/(^|f-), and 
fl = (^Hf")/(^if") m ^ ne 15- dimensional complex case. (This leads — multiplying by the 
corresponding i? 2 's, and || — to separability probabilities of ^ and J|, respectively.) 



2. Supplementary estimation of R\ constant 



In light of such imprecision, in our initial estimates, we undertook a supplementary 
analysis, in which, instead of examining each feasible 24-dimensional TF point for 2,001 
possible values of /i, with respect to separability or not, we simply used /i — 1. This, of 
course, allows us to significantly increase the number of points generated from the 64,000,000 
so far employed. 

We, thusly, generated 1,360,000,000 points, finding that we obtained a remarkably good 
fit to the important ratio Ri of the 24-dimensional measure (m sep ), at /i — 1, assigned to the 
separable two-qubit quaternionic density matrices to the (known) measure (m tot = 777 qq Q0 ) 
by setting Ri = (|f ) 2 ~ 0.114263 (our sample estimate of this quantity being the very close 
0.114262). This is exactly the square of the corresponding ratio || we had conjectured (based 
on extensive numerical and theoretical evidence) for the full (15-dimensional) complex two- 
qubit case in [17]. 
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3. Conjectured complex and quaternionic separability functions and probabilities 



Under this hypothesis on Ri for (3 = 4, we have the ensuing string of relationships 

sf u s M = (|) 2 (^ 3 -^>) 4 = (fd\^-^y= ( s "Um)\ (is) 

with (as already advanced in [17]), 

24/1 \ 2 6/ \ 2 

S^ex(jM) = ^(3-IM 2 )lM) =-((3-/i 2 )/ij . (14) 

Then, using our knowledge of the complementary ratio R 2 , given in (11), we obtain the 
desired exact result, 

r^Ht r, r, 72442944 2 6 -3 3 -7-53-113 n nrmn „„ 

= K,R, = = 52 . 17 . 19 . 23 . 7p » 0.0773765, (15) 

(the complex counterpart being as well as — in the framework of Andai [16] — that 

yHS = ovwit 26986 • 10~ 19 fie) 

sep/guat 24386773433626137413880000000 ~ ' K J 

B. Truncated quaternionic analysis ((3 = 3) 

For possible further insight into the HS two-qubit separability probability question, we 
undertook a parallel quasi-Monte Carlo (Tezuka-Faure) integration (setting fx = 1) for the 
truncated quaternionic case {(3 = 3), in which one of the four quaternionic parameters is 
set to zero. Although there was no corresponding formula for the HS total volume for this 
scenario given in [16], upon request, A. Andai kindly derived the result 

71-10 

Vt ™ nc = 384458588946432000 * 2 ' 43584 ' 10 ~ " 
In fact, Andai was able to derive one simple overall comprehensive formula, 

v" s = n^rx^ + 1) (is) 



yielding the total HS volumes for all n x n systems and Dyson indices (3. Let us, further, 
note that Andai obtains the result (17) as the product of three factors, V t ^ nc = i\\i\2^z-, 
where 

128tt 8 128 189tt 2 

7Ti = ; 7r 2 = ; 5T3 = . (19) 

105 893025 12696335643836416 v ; 
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Now, we will simply assume — in line with our basic Dyson-index ansatz, substantially 
supported in [17] and above — that the corresponding separability function is of the form 

« nc (/i)oc((3-/i 2 )/i)' 3 , = 3. (20) 

(Of course, one should ideally test this specific application of the ansatz too, perhaps using 
the same quasi-Monte Carlo method we have applied to the (3 = 4 instance above [Fig. 1].) 

We were somewhat perplexed, however, by the results of our quasi-Monte Carlo integra- 
tion procedure, conducted in the 18-dimensional space of off-diagonal entries of the truncated 
quaterionic density matrix p. Though, we anticipated (from our previous extensive numeri- 
cal experience here and elsewhere) that the estimate of the associated 18-dimensional volume 
would be, at least, within a few percentage points of = ^pf^ ~ 1-65793, our actual 
estimate was, in fact, close to 0.967 (1, thus, falling within the possible margin of error). 
Assuming the correctness of the analysis of Andai, which we have no other reason to doubt, 
the only possible explanations seemed to be that we had committed some programming 
error (which we were unable to discern) or that we had some conceptual misunderstanding 
regarding the analysis of truncated quaternions. (Let us note that we do convert the 4x4 
density matrix to 8 x 8 [complex] form [27, p. 495] [41, eq. (17)] [42, sec. II], while it 
appears that Andai does not directly employ such a transformation in his derivations.) 

In any case, we did devote considerable computing time to the (3 = 3 problem (generating 
1,180,000,000 18-dimensional Tezuka-Faure points), with the hope being that if we were in 
some way in error, the error would be an unbiased one, and that the all-important ratio of 
separable to total volume would be unaffected. 

Proceeding thusly, our best estimate (not making use of the Andai result (17) for the 
present) of the HS separability probability was 0.193006. One interesting possible candidate 
exact value is, then, ||| = 3 ^. 17 ~ 0.193062. (Note the presence of 128 in the numerators, 
also, of both factors tti and tt 2 , given in (19).) This would give us a counterpart [f3 = 3] 
value for the ratio R 2 of ^osreie ~ 0.278838. In [17], we had asserted that, in the other 
odd (3 = 1 case, the counterpart of R 2 was ^ril ~ 0.612315. (Multiplying this by ™§* 2 gave 
us the conjectured HS real two-qubit separability probability of ^.) 

So, let us say that although we believe we have successfully resolved — though still not 
having formal proofs — the two-qubit Hilbert-Schmidt separability probability question for 
the P = 2 and 4 (complex and quaternionic) cases, the odd (P = 1,3) cases, in particular 
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(3 — 3, appear at this point to be more problematical. 

C. Real and complex Qubit- Qutrit Hilbert- Schmidt Analyses 

For qubit-qutrit systems, we have previously reported [17, eq. (44)], following the lines 
of our (Bloore-parameterization-based) two-qubit analyses, that rather than the use of one 
ratio variable /i, in implementing the Peres-Horodecki positive-partial-transpose test for 
separability, it is necessary to employ two (corresponding specifically here to the case where 
the partial transpose is implemented by transposing the four 3x3 blocks of the 6x6 density 
matrix p in place) variables, already presented in (4). 

Once again, employing the Tezuka-Faure quasi-Monte Carlo methodology, we generated 
133,545 30-dimensional and 1,950,000 20-dimensional feasible points, corresponding now to 
the off-diagonal Bloore parameters of 6 x 6 complex and real density matrices, respectively. 
(Each analysis was run on a MacMini workstation for a number of weeks.) The much larger 
number of feasible Tezuka-Faure points generated in the real case was primarily due to our 
reparameterization in that case of the Bloore off-diagonal entries (essentially correlations) 
in terms of partial correlations [20-22] (cf. [46]). This allowed us to somewhat mitigate the 
computational "curse" of high dimensionality, in that each sampled point now corresponds 
to a density matrix and none (theoretically, at least) has to be discarded. (H. Joe has 
demonstated that it is possible to also implement this approach in the complex case, but the 
programming challenges for us were substantially greater, so we have not yet pursued such 
a course.) Of the 2,250,000 20-dimensional points, 38,622 were discarded because certain 
numerical difficulties (mainly convergence problems), arose in transforming to the partial 
correlations. The 133,545 feasible 30-dimensional points were drawn from (a much larger) 
430,000,000 ones. 

For each feasible sampled point we tested whether the associated (real or complex) 6x6 
density matrix was separable or not (that is, whether or not it passed the Peres-Horodecki 
test) for all possible pairs of v\ and u 2 ranging from to 1 in increments of j^-that is, 
101 2 = 10, 201 Peres-Horodecki positive-partial-transpose tests were performed for each 
feasible sampled Tezuka-Faure point. We present the two estimated bivariate separability 
functions in Fig. 3 and Fig. 4 (cf. [17, Figs. 3, 5]). 

In Fig. 5 we present a test of our Dyson-index HS separability-function ansatz by sub- 
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FIG. 3: Interpolated estimate over the unit square of the real qubit-qutrit separability function 
S^eai/qub-quti 171 ' 1 ' 2 )' based on 2,211,378 20-dimensional Tezuka-Faure points. For each of these 
points, 10,201 associated 6x6 density matrices, parameterized by v\ £ [0, 1] and vi £ [0, 1], were 
tested for separability. 




FIG. 4: Interpolated estimate over the unit square of the complex qubit-qutrit separability function 
^ 'complex /qub-qut^ 1 ' 1 ^ 1 ' 2 ^ based on 133,545 feasible 30-dimensional Tezuka-Faure points. For each 
point, 10,201 associated 6x6 density matrices, parameterized by v\ G [0, 1] and £ [0, 1], were 
tested for separability. 

tracting from Fig. 4 the square of the function in Fig. 3, which has been normalized so that 
its value at V\ = 1, = 1 equals that of the raw, unadjusted complex separability function. 
One should, of course, note the greatly-reduced z-axis scale from Fig. 4, indicating close 
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FIG. 5: The complex qubit-qutrit separability function shown in Fig. 4 minus the square of the 
real qubit-qutrit separability function (Fig. 4), the latter function normalized so that the value in 
the plot at v\ = 1, vi = 1 is 0. Note, importantly, the greatly- reduced z-axis scale vis-a-vis that of 
Fig. 4 

HS sep. funds. 



0.010 



0.005 




FIG. 6: The complex qubit-qutrit separability function (Fig. 4) and the normalized square of the 
real function (Fig. 3), holding v\ = i>i- By construction, the two curves are equal at v\ = i>2 = 1- 
The observed closeness of the two curves would be suggested by the HS Dyson-index ansatz 

adherence to the HS Dyson-index separability-function ansatz, which it has been a principal 
goal of this study to test. 

Now, in Fig. 6 we plot two very closely-fitting curves. One is the complex separability 
function holding v\ = i>2, and the other, the square of the real separability function also 
holding V\ = z/ 2 , but normalized to equal the first (complex) function at the point (1,1). 
This is also compelling evidence for the validity of the HS Dyson-index ansatz. 

Also, in [17], we indicated that it strongly appeared that though two ratio variables, v\ 
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FIG. 7: The estimated qubit-qutrit real separability function (Fig. 3) minus the candidate function 
(22), the latter function being scaled so the plotted value at (u±, v?) = (1, 1) is 0. 



and U2, given in (4), are ab initio necessary in the qubit-qutrit analysis, it seems that upon 
further analysis they coalesce into a product 

rj = Ul v 2 = , (21) 

P33P44 

and the separability function problem becomes actually simply univariate in nature, rather 
than bivariate. This aspect needs, of course, to be more closely evaluated in light of our 
new numerical results. In fact, one candidate HS separability function of such a univariate 
nature which can be seen to fit our estimated functions (Figs. 3 and 4) both very well (when 
appropriately normalized and/or squared) is (Fig. 7) 

S™i/ qub - q u t M oc 1 - (1 - = 1 - (1 - 77)1 = h -B v {l, h -). (22) 

In Fig. 7 we show the fit of this function to the estimated qubit-qutrit real separability 
function (Fig. 3). In Fig. 8 we show the sum-of-squares of the fit of the one-parameter 
family of functions 1 — (1 — V\V-ip to the normalized estimated real qubit-qutrit separability 
function (Fig. 3). (For 7 = | we obtain (22). We observe that the minimum of the curve is 
somewhat in the neighborhood of 7 = |.) 

In Fig. 9 we show the sum-of-squares of the fit of the two-parameter family of functions 
1 — (1 — (yxV'ifY 1 to the normalized estimated real qubit-qutrit separability function (Fig. 3). 
(For 7 = |, 9 — 1 we obtain (22). We observe that the minimum of the curve is somewhat 
in the neighborhood of 7 = |, 9 = 1.) 
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FIG. 8: The sum-of-squares (SS) of the fit of the one-parameter family of functions 1 — (1 — v\V2p 
to the normalized estimated real qubit-qutrit separability function (Fig. 3). For 7 = 2.5, we obtain 
the candidate separability function (22) 




FIG. 9: The sum-of-squares (SS) of the fit of the two-parameter family of functions 1 — (1 — (z^i^) 9 ) 7 
to the normalized estimated real qubit-qutrit separability function (Fig. 3). For 7 = 2.5, 6 = 1, we 
obtain the candidate separability function (22) 



In the framework of Andai [16], the total volume of the 25-dimensional convex set of real 
6x6 density matrices can be represented as the product 



V, 



HS 



8192tt 6 



25tt 3 



7T 



qub-qut/reai 2 53125 1399771004732964864 1730063650258944000 



1.72301 ■ 10 



-ii 



(23) 



and the total volume of the 35-dimensional convex set of complex 6x6 density matrices as 
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the product 

TT 15 I 



qub-qut/ complex 86400 000 3460550346681745424512204800 V ; 

71-15 

9.58494 • 10" 29 . 



298991549953302804677854494720000000 
In both of these products, the first (20- or 30-dimensional) factor serves as the denominator 
of the ratio R\ and the second (5-dimensional) factor, as the denominator of the ratio R 2 . 
(The numerator of R\ is the 20- or 30-dimensional mass assigned to the separable density 
matrices, and the numerator of R2 is the integral of the product of the corresponding suitably 
normalized separability function and Bloore jacobian over the 5-dimensional unit simplex.) 

Under the assumption of the correctness of (22), we find that the qubit-qutrit counterparts 
of the important constants R2 employed above in the two-qubit case ((11) and immediately 
below there) are 



«V,_„,(/3 = 1) = 1 - ,,„„„„ ,, r ^ ~ 0.724715, (25) 
«W*0> = 2 > = ~™^7Z:Z r^"'"^' - 0.692789, (26) 



4194304 
4849845tt 
44632342463 + 68578836480 log(2) 
4190140110 

^V 6 _ 9Ut 03 = 4)= (27) 

192210846322598002116984324520591 - 277301145703236210250598232096768 log(2) 

501570554133080277487570824 

« 0.675902, 

for the real, complex and quaternionic cases, respectively. Again, we note a monotonic 
decrease as the Dyson index (5 increases. Also, for the truncated quaternionic {(5 = 3) case, 
we found 

^V 6 _ 9Ut (/3 = 3)= (28) 

_ 967504709 18446744073709551616(-67294453713397888 + 5638997741091tt) 
552123 71729672378917671400466262753675tt 2 

« 0.681261. 

Our sample estimates of the complementary R 1 constants were 0.226468, in the real case, 
and 0.047679 in the complex case. (We note that these estimates should be independent of 
the choice (correctness) of separability functions (22) and, in the complex case, the square 
of (22).) Forming the products R±R 2 based on these estimates, we obtain an estimated real 
separability probability of 0.164125 and complex separability probability of 0.0330446. 
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Our numerical analyses have been concerned only with pairs of values of v\ and u 2 lying 
within the (bounded) unit square [0, 1] x [0,1]. We have made the implicit assumption that 
the substance of our analyses would not be altered/biased if we were able to incorporate 
all possible pairs lying within the unbounded quadrant [0, oo] x [0, oo]. (Clearly, for points 
v\ > l,i/2 > 1) simply by symmetry considerations, we immediately expect the separability 
function to be proportional to 1 — (1 — For pairs of points (ui, z/ 2 ) for which one member 
is greater than 1 and the other less than one, we expect the form the separability function 
takes to depend on whether or not rj — V]V 2 > !■) 

1. Relations to previous qubit-qutrit analyses [11, 17] 

In [9], we had undertaken a large-scale numerical (again, Tezuka-Faure quasi-Monte Carlo 
integration) analysis of the separable volumes of the 35-dimensional convex set of complex 
qubit-qutrit systems, endowed with the Hilbert-Schmidt, as well as a number of monotone 
(including the Bures) metrics. The estimate we obtained there of the HS separability prob- 
ability was 0.0266891. As we pointed out in our subsequent study [17, p. 14305], this is 
remarkably close to ^ « 0.0266889. 

In the subsequent study [17], which was chiefly devoted to the case of two-qubit systems, 
we had included supplementary analyses of the real and complex qubit-qutrit systems. But 
there we had only employed-in the interest of alacrity-Monte Carlo (random number), rather 
than ("lower-discrepancy") quasi-Monte Carlo methods. So, the results presented in this 
paper, we believe should be more accurate and informative. (Also, rather than sampling 
grids for v-i,v 2 of size 101 x 101, we had employed grids of sizes 50 x 50 in the real case, 
and 20 x 20 in the complex case.) In [17, sec. 10], we had put forth the tentative hypothesis 
that the corresponding real qubit-qutrit separability function was simply proportional to 
^/fj = ^Jv\v 2 (having a somewhat similar profile to our present candidate (22), but definitely 
providing an inferior fit to our numerical results here). 

D. Bures two-qubit separability functions and probabilities 

Now, in our present study, we shall somewhat parallel the sequential approach of 
Zyczkowski and Sommers in that they, first, computed the total volume of (separable and 
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nonseparable) n x n density matrices in terms of the (flat or Euclidean) Hilbert-Schmidt 
metric [13] [47, sees. 9.6-9.6, 14.3], and then, using the fundamentally important Bures 
(minimal monotone) metric [47, sec. 14.4] [14]. (In particular, they employed the Laguerre 
ensemble of random matrix theory [15] in both sets of computations (cf. [16]). The Bu- 
res and HS metrics were compared by Hall [48], who concluded that the Bures induced 
the "minimal- knowledge ensemble" (cf. [49]), also noting that in the single-qubit case, the 
Bures metric "may be recognized as the spatial part of the Robertson- Walker metric in 
general relativity".) That is, we will seek now to extend the form of analysis applied in the 
Hilbert-Schmidt context in [17] to the Bures setting. 



1. Review of earlier parallel Hilbert-Schmidt findings 

To begin, let us review the most elementary findings reported in [17, sec. II.A.l]. The 
simplest (four-parameter) scenario studied there posits a 4 x 4 density matrix p with fully 
general diagonal entries (pu, P22, P33, P44 = 1 — P11 — P22 — P33) and only one pair of real 
off-diagonal non-zero entries, P23 = P32- The HS separability function for that scenario was 
found to take the form [17, eq. (20)], 

S[(£)]M=< ' ' • (29) 




where we primarily employ the variable p = y ; rather than v — p 2 , as in [11, 17]. 

Allowing the 23- and 32-entries to be complex conjugates of one another, we further found 
for the corresponding separability function [17, eq. (22)] — where the wide tilde over an i,j 
pair will throughout indicate a complex entry (described by two parameters)— 



^M-ty^M)'-*"'" (30) 



Trp 2 < p < 1 

7T p > 1 



Further, permitting the 23- and 32-entries to be quaternionic conjugates of one another 
[27, 50], the corresponding separability function [17, eq. (24)] — where the wide hat over an 
i,j pair will throughout indicate a quaternionic entry (described by four parameters) — took 
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the form 



^ < fj, < 1 



2 



/i > 1 



So, the real (29), complex (30), and quaternionic (31) HS separability functions accord 
perfectly with the Dyson index sequence f3 = 1,2,4 of random matrix theory [33]. "The 
value of j3 is given by the number of independent degrees of freedom per matrix element and 
is determined by the antiunitary symmetries ... It is a concept that originated in Random 
Matrix Theory and is important for the Cartan classification of symmetric spaces" [51, p. 
480]. The Dyson index corresponds to the "multiplicity of ordinary roots" , in the terminology 
of symmetric spaces [52, Table 2]. However, we remain unaware of any specific line of 
argument using random matrix theory [15] that can be used to formally confirm the HS 
separability function Dyson-index-sequence phenomena we have noted above and observed 
in [17] . (The basic difficulty/novelty appears to be that the separability aspect of the problem 
introduces a totally new set of complicated constraints — quartic (biquadratic) in /j, [17, eq. 
(5)] [11, eq. (7)] — that the multivariate integration must respect [17, sec. I.C].) 

As a new exercise here, unreported in [17], we found that setting any single one of the 
four components of the quaternionic entry, x 2 3 + W23 + jjW23 + ki> 2 3, in the scenario just 
described, to zero, yields the (truncated quaternionic) separability function, 



4:TTU 



3 



< n < 1 



'In 



consistent, at least, in terms of the exponent of /i, with the Dyson-index pattern previously 
observed. 

Continuing the analysis in [17], we computed the integrals 

POO 

V HS = / <*HS / \ jHS ( \ 1 /oo\ 

" sep I scenario I scenarioKr 1 I^J scenario Kr 1 j^f^i V 00 / 

JO 

of the products of these separability functions with the corresponding (univariate) marginal 
jacobian functions (which are obtained by integration over diagonal parameters only and 
not any of the off-diagonal rc^'s and y^'s) for the reparameterization of p using the Bloore 
variables [17, eq. (17)]. This yielded the HS scenario-specific separable volumes V s ^p/ scenario - 
The ratios of such separable volumes to the HS total volumes 

POO 

■io / JscenarMdfi, (34) 
JO 



f oo 

yHS = HS 

tot I scenario ^scenario 
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where c^ nario is a scenario-specific constant, gave us in [17] (invariably, it seems, exact) 
separability probabilities. (For the three scenarios listed above, these probabilities were, 
respectively, ||, | and ^.) 

Based on the numerous scenario-specific analyses in [17], we are led to believe that the 
real, complex and quaternionic separability functions conform to the Dyson-index pattern for 
general scenarios, when the Hilbert-Schmidt measure has been employed. This apparent ad- 
herence was of central importance in arriving at the assertions in [17, sees. IX. A. 1 and IX. A. 2] 
that the HS separability probabilities of generic (9-dimensional] real) and (15-dimensional) 
complex two-qubit states are ^ and J|, respectively. There we had posited — using mutually 
supporting numerical and theoretical arguments — that [17, eq. (102)] 

«5&)oc^(3-/,>, (35) 

and, further pursuing our basic Dyson-index ansatz (fitting our numerical simulation ex- 
tremely well [17, Fig. 4]), that {S^{p)) 2 oc S^ plex {p). (Also, in the first part of the 
analyses above, we presented numerical evidence that (<S^(//)) 4 oc S^f at {p), and made this 
relation more precise (13).) 



2. Four-parameter density-matrix scenarios 

Now, employing formulas (13) and (14) of Dittmann [53] for the Bures metric — which 
avoid the possibly problematical need for diagonalization of p — we were able to find the Bures 
volume elements for the same three basic (one pair of free off-diagonal entries) scenarios. 
We obtained for the real case, 



;»,.,,, y/Pll V 1 ~ Pll _ P22 

dV [(2>3)] = — — -==== dp u dp 22 dx 23 dp, (36) 

VI - x 2 3 (p 22 /i 2 + p u ) vVp| 2 + (1 - Pu) Pn 

for the complex case, 



7t / Bures PllP22 (Pll + P22 ~ 1) A A A A A (Q>7\ 

dv \r^\ = ,/i 2 ^r, — ; w — i ^Y\ dpildp22dX2 * dy2zd ^ 

V 1 - Vh - X h (P22P 2 + Pll) (-P?1 + Pll + P 2 P 2 22 ) 

and for the quaternionic case, 

dV [j2§j] = ^ d P^ d P^ dx ^dy 2 3du 23 dv 23 dp, (38) 

where 

A = -P\\P%2 (Pll + P22 - I) 2 , 
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and 

B = 4^1 - w| 3 - v 2 3 -x 2 23 -y 2 3 (p 22 p 2 + Pn) (-Pn + Pn + p' l p 2 22 f ■ 

In analyzing the quaternionic case, we transformed — using standard procedures [27, p. 
495] [41, eq. (17)] — the corresponding 4x4 density matrix into an 8 x 8 density matrix with 
(only) complex entries. To this, we found it most convenient to apply — since its eigenvalues 
and eigenvectors could be explicitly computed — the basic formula of Hiibner [54] [55, p. 
2664] for the Bures metric. 

Integrating these three volume elements over all the (four, five or seven) variables, while 
enforcing the nonnegative definiteness requirement for p, we derived the Bures total (sep- 
arable and nonseparable) volumes for the three scenarios — V t ot/[(23)] = ~ 0.822467, 
V Burej_ = £ w 0.484473, and V Bw 2~~ = « 0.126835. 



4 

VI 

tot/[(2,3)] 64 ' ' tot/[(2,3)} ~ 

We note importantly that the Bures volume elements ((36), (37), (38)), in these three 



cases, can be factored into products of functions of the off-diagonal Bloore variables, 
M 23? ^23) ^23 an d U23-, and functions of the diagonal variables, pu,p 22 and p. Now, we 
will integrate (one may transform to polar and spherical coordinates, as appropriate) just 

the factors , 1 = , , = = and , 3 = = = — involving the off-diagonal vari- 

able(s) over those variables. In doing this, we will further enforce (using the recently- 
incorporated integration-over-implicitly-defined-regions feature of Mathematica) the Peres- 
Horodecki positive-partial-transpose-criterion [29, 30, 56], expressible as 



in the real case, 



in the complex case, and 



p 2 -x 2 23 >0 (39) 



^ - 4s ~ 2/23 > 0, (40) 



- ^23 - 2/23 - M 23 - *4 > °> ( 41 ) 

in the quaternionic case. (None of the individual diagonal pa's appears explicitly in these 
constraints, due to an attractive property of the Bloore [correlation coefficient /off-diagonal 
scaling] parameterization. Replacing p 2 in these three constraints by simply unity, we obtain 
the non-negative definiteness constraints on p itself, which we also obviously must enforce.) 
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Performing the indicated three integrations, we obtain the Bures separability functions, 



7T /i > 1 

2sin~ 1 ( J u) < ji < 1 



$MfM = < ' " " ( 42 ) 




u > 1 

Sfj^U) = { , (43) 

[(2,3)] 1 rt ^ /- -\ V ^ 

- a/I - At 2 ) < // < l 



and 



47T 

Sp^U) = \ 3 ' . (44) 

1 |7r 2 (- v ^7?/i 2 -2 v ^7? + 2) 0<fi<l 

Then, utilizing these three separability functions — that is, integrating the products of the 
functions and the corresponding remaining diagonal-variable factors in the Bures volume 
elements ((36), (37)), ((38)) over the fi,pn and P22 variables — we obtain separable volumes 
of V£^ 3)] = 0.3658435525 and 

yBures = y^'-es _ 2/_ 2C - + n \ = — 7r 2 (4C - 6 + tt) « 0.124211 (45) 

«ep/[(2,3)] tot/[(2,3)] 32 V 7 64 V 7 V 7 

and consequent separability probabilities, respectively, of 0.4448124200 and (our only exact 
Bures separability probability result in this study (cf. [6])), 

pBures = 4C - 6 + 7T _ Q.256384, ( 46 ) 

«ep/[(2,3)] 7T ' K J 

where C ~ 0.915966 is Catalan's constant (cf. [57]). (This constant appears commonly in 
estimates of combinatorial functions and in certain classes of sums and definite integrals [58, 
sec. 1.7]. The ratio ^-as well as having an interesting series expansion [58, p. 54]-occurs in 
exact solutions to the dimer problem of statistical mechanics [58, p. 54] [59, 60]). 

Further, for the quaternionic case, V Bur S^ « 0.012954754466, and p Bur f^ « 

' •cp/[(2,3)] «ep/[(2,3)] 

0.10213883862. (The corresponding HS separability probability was also of the same rela- 
tively small magnitude, that is, ^ [17, sec. II.A.3]. We have computed the various Bures 
separable volumes and probabilities to high numerical accuracy, hoping that such accuracy 
may be useful in searches for possible further exact formulas for them.) 

So, the normalized — to equal 1 at /i — 1 — forms of these three separability functions 
are — ^ i(2,3)] — an( j — — _ j n pjg^ we pi Q t — motivated by the appearance of 

the Dyson indices in the analyses of [17] — the fourth power of the first (real) of these three 
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At > 1 




FIG. 10: Joint plot of the normalized Bures quaternionic separability function — — , the 
square of the normalized Bures complex separability function [(2 '|^ — , and the fourth power of the 
normalized Bures real separability function [(2|3)1 . 



The order of dominance of the three curves 



is the same as the order in which they have been mentioned. 

normalized functions together with the square of the second (complex) function and the 
(untransformed) third (quaternionic) function itself. We find a very close, 



Srfo U ^ S Gu)\4 



1(2,3)] 



7T 



£ Bures ( iA 

" 2it ' 



- [(2,3)] ^' 
47T 2 



), 



(47) 



but now not exact fit, as we did find in [17] for their (also normalized) Hilbert-Schmidt 
counterparts ^f™, -^f^- and ^ ((29), (30), (31)). 

As an additional exercise (cf. (32)), we have computed the Bures separability function 
in the (truncated quaternionic) case that a single one of the four components of the (2,3)- 
quaternionic entry is set to zero. Then, we have (falling into the same tight cluster in Fig. 10, 
when the ^-power of its normalized form is plotted) 



cBures 



§tt 2 (4- V21og(3 + 2V2)) n>l 
\tx Uy/1 - H 2 - sin-^/i)) (y/2 log (3 + 2^2 - 4)) < y, < 1 



(48) 



We have been able, further, using the formulas of Dittmann [53] , to compute the Bures 
volume elements for the corresponding (five-dimensional) real and (seven-dimensional) com- 
plex scenarios, in which both the {2,3} and {1,2} entries are allowed to freely vary. But 
these volume elements do not appear, now, to fully factorize into products of functions (as 
is the case for (36) and (37)) involving just pn, p 2 2, P and just the off-diagonal variables x^'s 
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and i/i/s. The requisite integrations are, then, more problematical and it seemed impossible 
to obtain an explicit univariate separability function of /i. 

For instance, in this regard, we have for the indicated five-dimensional real scenario that 



1 / A 

^^(2,3)] = 1\I BC ^ D + E ^ d Pll dp22dx 12 dx 23 dfi, (49) 

where 

A = -Pup1 2 (pn + P22 - 1) ( {p 2 - 1) P22 + 1) , (50) 
B = (P22P 2 + Pn) 2 ,C = x\ 2 + xj 3 - 1, (51) 
D = (pn + P22) (x\ 2 p 22 (P22P 2 + Pn) 2 ~ ((p 2 ~ 1) P22 + 1) {-Pu + Pn + P 2 pl 2 )) (52) 
and 

E = ~x\zP22 (pn + P22 ~ 1) (-pn + pn + P 2 p\ 2 ) ■ (53) 
So, no desired factorization is apparent. 

3. Five-parameter density-matrix scenarios 

However, the computational situation greatly improves if we let the (1,4) and (2,3)- 
entries be the two free ones. (These entries are the specific ones that are interchanged 
under the operation of partial transposition, so there is a greater evident symmetry in such 
a scenario.) Then, we found that the three Bures volume elements all do factorize into 
products of functions of off-diagonal entries and functions of diagonal entries. We have 



dv ^l)kn - - {x 2 u _ 1} {X 2 3 _ 1} {p ] 2 + ps3 _ 1} (p22 + p3s) d Pll d P22 d P33 dx 14 dx 23 , 

(54) 

where simply for succinctness, we now show the volume elements before replacing the p 33 
variable by ji. (We note that the expression for dV^ u ^ 23 ^ is independent of pn.) For the 
corresponding complex scenario, 

dV ixim] = ^/I d Pll d P22 d P33 dr l4 dr 23 d9 u d9 23 , (55) 

where 

F = -rl 4 rl 3PuP 22p33 (pn + P22 + P33 ~ 1) , (56) 
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and 



G = (r 2 4 - 1) (r 2 23 - 1) (p 22 + pas - I) 2 (P22 + Pss) 2 , 



(57) 



and we have now further shifted to polar coordinates, Xij + iy^ - , ;)i . ; 



i sin 6ij) 



T ij [COS ( 

For the quaternionic scenario, we have (using two sets of hyperspherical coordinates 

(r 14, &u > #14) #14 ) an d (?" 23, #23 > ^23^ ^))) 



j-tr Bures 
[(1,4),(2,3)] 



1 



^dp^dp^df^dr^T^dB^dB^dB^ (58) 



where 



F = sin 2 (6fi>) sin (flff) sin 2 sin (*g>) r^fV^ (- Pll - P22 - P33 + l)* 2 

(59) 

and 

^ ^ ' (60) 



G 



1 — r 



Az (P22 + P33 - 1) (P22 + p3sY 



The total Bures volume for the first (real) of these three scenarios is V^j*^ ^ 2 
0.484473, for the second (complex) scenario, V Bur ^- = 77^ 

tot/[(l,4),(2,3)] 192 



3)] 64 



0.507339, and for the 



third (quaternionic), y Bur fS~^ = ^ 

VH y ' tot/[(l,4),(2,3)] 245760 



0.0039119. 



In the two corresponding Hilbert-Schmidt (real and complex) analyses we have previously 
reported, we had the results [17, eq. (28)], 



5 [fl 5 4),(2,3)](p) 



and [17, eq. (34)] 



<i HS (iA 



4p < fj, < 1 
p > 1 



7T V < P < 1 



(61) 



(62) 



thus, exhibiting the indicated exact (Dyson-index sequence) proportionality relation. We 
now found, for the two Bures analogs, that 



5 [fMM2,3)](/ i ) - < 



7T 2 1 
27T CSC~ 1 (p) /i > 1 

27rsin~ 1 (/i) < p < 1 



(63) 
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H = l 




i67r 2 



167T 2 (1- V / W^ 



/i > 1 



< /i < 1 



(64) 



and, further still, for the quaternionic scenario, 




j7r 4 (- v ^7^ /U 2 -2 v ^^ + 2) 0</i<l 




and quaternionic separability probabilities of 0.3042243652, 0.19102778 and 0.120436049. 

We see that for values of /i e [0, 1], the normalized forms of these three Bures separability 
functions are identical to the three obtained above ((42), (43), (44)) for the corresponding 
single- nonzero-entry scenarios. While those earlier functions were all constant for /i > 1, we 
now have symmetrical behavior about n = 1 in the form, Sf^ZM = ^ceZrioij,)- 

In Fig. 11, we show the analogous plot to Fig. 10, using the normalized (to equal 1 at 
ji — 1) forms of the three additional Bures separability functions ((63), (64), (65)). We 
again, of course, observe a very close fit to the type of proportionality relations exactly 
observed in the Hilbert-Schmidt case ((61), (62)). 

We were, further, able to compute the Bures volume element for the £/iree-nonzero-entries 
complex scenario [(1, 2), (1, 4), (2, 3)], but it was considerably more complicated in form than 
those reported above, so no additional analytical progress seemed possible. 

4- Additional analyses 

Regarding the possible computation of Bures separability functions for the totality of 
9-dimensional real and 15-dimensional complex two-qubit states, we have found, prelimi- 
narily, that the corresponding metric tensors (using the Bloore parameterization [sec. I A]) 
decompose into 3x3 and 6x6, and 3x3 and 12 x 12 blocks, respectively. The 3x3 blocks 
themselves are identical in the two cases, and of precisely the (simple diagonal) form (if we 
employ hyperspherical coordinates) that Akhtarshenas found for the Bures metric using the 
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FIG. 11: Joint plot of the normalized Bures quaternionic separability function — [(1,4) ^ 2|3)1 — , the 
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square of the normalized Bures complex separability function 1(1 ' 4 jg^2 3)1 — > an d the fourth power of 

^Bures ((-l') 

the normalized Bures real separability function EA(|ig)] — _ Over the interval \i £ [0, 1], the three 
functions are identical — with the same order of dominance — to those in Fig. 10. 

coset parameterization [61, eq. (23)]. These 3x3 blocks, thus, depend only upon the diag- 
onal entries (while in [61], the dependence, quite differently, was upon the eigenvalues). It 
appears, though, that the determinants — for which we presently lack succinct formulas — of 
the complementary 6x6 and 12 x 12 blocks, do depend upon all, diagonal and non-diagonal, 
parameters, rendering further analytical progress, along the lines pursued with substantial 
success for the Hilbert-Schmidt metric, for these scenarios rather problematical. 

5. Discussion 

The close proximity observed above between certain two-qubit separability results for 
the Hilbert-Schmidt and Bures metrics is perhaps somewhat similar in nature/explanation 
to a form of high similarity also observed in our previous analysis [9]. There, large scale 
numerical (quasi-Monte Carlo) analyses strongly suggested that the ratio of Hilbert-Schmidt 
separability probabilities of generic (rank-6) qubit-qutrit states (6x6 density matrices) to 
the separability probabilities of generically minimally degenerate (boundary/rank-5) qubit- 
qutrit states was equal to 2. (This has since been formally confirmed and generalized — in 
terms of positive-partial-transpose-ratios — to arbitrary bipartite systems by Szarek, Bengts- 
son and Zyczkowski in [62]. They found that the set of positive-partial-transpose states is 
"pyramid decomposable" and, hence, is a body of constant height.) Parallel numerical ratio 
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estimates also obtained in [9] based on the Bures (and a number of other monotone) metrics 
were also surprisingly close to 2, as well (1.94334 in the Bures case [9, Table IX]). However, 
no exact value for the Bures qubit-qutrit ratio has yet been established, and our separability 
function results above, might be taken to suggest that the actual Bures ratio is not exactly 
equal to 2, but only quite close to it. (Possibly, in these regards, the Bures metric might 
profitably be considered as some perturbation of the flat Euclidean metric (cf. [63]).) 

Further study of the forms the Bures separability functions take for qubit-qubit and 
qubit-qutrit scenarios is, of course, possible, with the hope that one can gain as much insight 
into the nature of Bures separability probabilities as has been obtained by examining the 
analogous Hilbert-Schmidt separability functions [17]. (In [8], we had formulated, based 
on extensive numerical evidence, conjectures — involving the silver mean, y/2 - 1— for the 
Bures [and other monotone metric] separability probabilities of the 15-dimensional convex 
set of [complex] qubit-qubit states, which we would further aspire to test. One may also 
consider the use of monotone metrics other than the minimal Bures one [16] — such as the 
Kubo-Mori and Wigner-Yanase.) The analytical challenges to further progress, however, in 
light of the apparent non-factorizability of volume elements into diagonal and off-diagonal 
terms, appear quite formidable. 

III. EULER- ANGLE-PARAMETERIZATION SEPARABILITY FUNCTIONS 

In the previous section (sec. II D), we found that the Bloore parameterization (sec. I A), 
markedly useful in determining separable volumes based on the (non-monotone) Hilbert- 
Schmidt metric, is less immediately fruitful when the Bures (and presumably other mono- 
tone) metrics is employed. In light of this development, it appeared to be of interest to see if 
some other parameterizations might prove amenable to such type of "separability function" 
analyses. In particular, we will examine here the use for such purposes-as earlier suggested 
in [34]-of the SU (4) Euler-angle parameterization of the 15-dimensional complex set of 4 x 4 
density matrices developed by Tilma, Byrd and Sudarshan [31]. (We will closely follow the 
notation and terminology of [31]. In [34], we simply attempted to fit symmetric polynomials 
[64] to yield previously-conjectured separable volumes, and did not initiate any independent 
quasi-Monte Carlo estimation procedures, as we will immediately below.) 
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A. Complex two-qubit case 

The fifteen parameters, then, employed will be twelve Euler angles (ctj's) and three in- 
dependent eigenvalues (Ai, A 2 , A 3 ). The total (separable and nonseparable) volume is simply 
(for all metrics of interest) the product of integrals over these two sets of parameters [13, 14]. 
Now, to study the separable-volume question, in complete analogy to our methodology 
above, we would like to integrate over the twelve Euler angles (rather than the off-diagonal 
Bloore parameters, as before), while enforcing the positivity of the partial transpose, required 
for separability. We were, fortunately, able to perform such enforcement in the Bloore- 
parameterization case employing only a single diagonally-related parameter (/z), given in 
(3). Such a reduction in the number of relevant parameters, however, does not seem possi- 
ble in the Euler-angle parameterization case. So, the analogue of the separability function 
we will obtain here will be a trivariate function (of the three eigenvalues). Hopefully, we will 
be able to determine an exact functional form for it, and then utilize it in further integra- 
tions to obtain separable volumes, in terms of both monotone and non-monotone metrics. 
(Also, the question of whether counterpart Euler-angle separability functions in the real and 
quaternionic cases adhere to some form of Dyson-index-sequence behavior certainly merits 
attention.) 

B. Trivariate separability function for volume computation 

Now, we made use of a sequence of 1,900,000 12-dimensional Tezuka-Faure points (twelve 
being the number of Euler angles over which we will integrate). For each such point, we let 
the associated three (free) eigenvalues each take on all possible values from ^ to 1, in steps 
of 4j. Of course, the possible triads of free eigenvalues is constrained by the requirement 
that they not all sum to more than 1. There were 9,880 such possible triads. For each such 
triple-holding the Euler angles constant-we evaluated whether the associated 4x4 density 
matrix was separable or not. 

We, then, interpolated the results to obtain functions defined over the three-dimensional 
hypercube [0, l] 3 . In Fig. 12, we show a two-dimensional marginal section (over Ai,A 2 , 
say) of this function (obtained by summing over the values of A 3 ) of the estimated three- 
dimensional separability function. We know from the work of Pittenger and Rubin [65, Cor. 
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FIG. 12: Two-dimensional marginal section of the estimated three-dimensional separability func- 
tion based on the Euler- angle parameterization for the 15-dimensional convex set of complex 4x4 
density matrices. Note the mesa/plateau shape, indicative of the fully separable neighborhood of 
the fully-mixed state 

4.2], for example, that, for the specific case of two-qubits, any density matrix all of the 
four eigenvalues of which are greater than ^ w 0.2333 must be separable. Therefore, there 
certainly does exist a neighborhood of the fully-mixed state (having its four eigenvalues equal 
to j) that is composed of only separable states. This is reflected in the plateau present in 
Fig.. 12. 

In Fig. 13 we, additionally display a one-dimensional marginal section (over Ai), obtained 
by summing over both Ai and A2, of the estimated three-dimensional separability function. 
(The curve now appears unimodal rather than flat at its maximum,) 

Employing the trivariate separability function obtained by interpolation from the data 
we have generated here, we were able to obtain an estimate of 0.242021 for the HS sep- 
arability probability. (From our extensive Bloore-parameterization analyses, as previously 
noted we believe its true value is J| ~ 0.242424.) We were essentially just as readily able to 
obtain an estimate of the Bures (minimal monotone) separability probability (or any of the 



33 



marc in a] &cp. ftmcL 




— , i , ; — i — i , 1 i — , , 1 , , , 1 — — , 1_ x 

0.4 0.6 O.S ]JQ 

FIG. 13: One-dimensional marginal section of the estimated three-dimensional separability function 
based on the Euler-angle parameterization for the 15-dimensional convex set of complex 4x4 density 
matrices 

other monotone metrics-Kubo-Mori, Wigner-Yanse,. . . [66], it appears) of 0.0734223, while 
in [8], this had been conjectured to equal 1680 (^ +v ^) ~ 0.733389. (In the HS and Bures 
computations reported here and in the next section, we perform numerical integrations over 
the simplices of eigenvalues, in which the integrands are the products of our interpolated 
separability functions and the appropriate scenario-specific volume elements indicated in the 
twin Sommers-Zyczkowski 2003 papers [13, 14].) 

Of course, now the research agenda should turn to the issue of finding a possibly exact 
formula (undoubtedly symmetric in the three eigenvalues [34, sees. V and VI]) for this 
three-dimensional Euler-angle-based separability function, and for other qubit-qubit and 
qubit-qutrit scenarios. 

In Fig. 14 we show, based on the 9,880 points sampled for each 12-dimensional TF-point, 
the estimated value of the separability function for that point paired with the Euclidean 
distance of the vector of eigenvalues (Ai, A2, A3, 1 — Ai — A2 — A3) for that point from the 
vector of eigenvalues (7, 7, 7, 7), corresponding to the fully mixed state. 
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triv. sep. funct 




FIG. 14: The trivariate Euler-angle complex qubit-qubit separability function paired with the Eu- 
clidean distance of the corresponding vector of eigenvalues to the vector {\, \ , \ , \ ), corresponding 
to the fully mixed state 

C. Bivariate separability function for area computation 

Now we repeat the procedures described immediately before, except for the a priori 
setting of one of the three free eigenvalues to zero, so the associated density matrices must 
lie on the 14-dimensional boundary of the 15-dimensional convex set of two-qubit complex 
states. (The analysis was conducted independently of that pertaining to the volume, and 
now we were able to employ a much larger number-23,500,000-of TF-points.) The resulting 
separability function is now bivariate, lending itself immediately to graphic display. In 
Fig. 15 we show this function, and in Fig. 16, its one-dimensional section over Ai. 

Employing the bivariate separability function (Fig. 15) obtained by interpolation from the 
data (23,500,000 TF-points) we generated, we were able to obtain an estimate of 0.12119 for 
the HS separability probability, which, from our complementary Bloore analyses, together 
with the ("one-half") Theorem 2 of [62], we believe to be exactly ^ « 0.121212. (The 
proximity of our estimate to this value clearly serves to further fortify our conjecture that 
the HS separability probability of generic complex two-qubit states is J|.) Additionally, our 
estimate of the associated Bures separability probability was 0.0396214, approximately one- 
half that of the corresponding probability for the non-degenerate complex two-qubit states 
[9]. 

In Fig. 17 we show, based on 780 points sampled for each 12-dimensional TF-point, the 
estimated value of the bivariate separability function for that point paired with the Euclidean 
distance of the vector of eigenvalues for that point from the vector of eigenvalues (|, |, |, 0) 
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FIG. 15: Bivariate Euler-angle separability function for the 14-dimensional boundary hyperarea of 
the 15-dimensional convex set of complex 4x4 density matrices 

marginal uep. fund. 




FIG. 16: One-dimensional marginal section of the estimated two-dimensional Euler-angle-based 
separability function for the 14-dimensional hyperarea of the 15-dimensional convex set of complex 
4x4 density matrices 
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FIG. 17: Paired values of the estimated bivariate Euler-angle separability function (Fig. 15) and the 
Euclidean distance of the associated vector of eigenvalues to the vector (|, |, |,0) corresponding 
to the most fully mixed boundary state 



corresponding to the most fully mixed boundary state. 



R(p) = — . (66) 



D. Participation ratios 

The participation ratio of a state p is denned as [1, eq. (17)] [47, eq. (15.61)] [67, 68] 

1 

trp 2 

For R(p) > 3, a two-qubit state must be separable. For convenience, we will also employ 
the variable 

s w " " W)>' (67) 

which varies over the interval [0,1] for states outside the separable ball. In Fig. 18 we show a 
plot-having set one of the four eigenvalues to zero-of the sixth-power of this ratio. We note 
a close similarity in its shape to the estimated bivariate separability function displayed in 
Fig. 15. Although we can not similarly visually display the trivariate separability function 
we have found that it is closely fit-outside the separable ball (R(p) > 3) [47, Fig. 15.7]- by 
the fourth-power of the participation ratio. In Fig. 19, we show the trivariate separability 
function vs. R(p) 4 . (For R(p) 4 > 81, only separable states are encountered). In Fig. 20, we 
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FIG. 18: The sixth-power, R(p) G , of the participation ratio for minimally degenerate two-qubit 
states with one eigenvalue equal to zero. Note the similarity to Fig. 15 

triv. scp . fiance 



FIG. 19: The (Euler-angle) trivariate separability function for generic (15-dimensional) complex 
two-qubit states plotted against the fourth-power of the participation ratio for the 9,880 points 
sampled 

show the comparable bivariate separability function vs. R(p) 6 - 

As an exercise (not directly tied to our quasi-Monte Carlo computations), we assumed 
that the trivariate Euler-angle separability functions are all proportional (outside the sepa- 
rable ball, R(p) > 3) to some powers of R(p), for each of the generic real, complex, truncated 
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FIG. 20: The (Euler-angle) bivariate separability function for generic (14-dimensional) minimally 
degenerate complex two-qubit states plotted against the sixth-power of the participation ratio for 
the 780 points sampled 

quaternionic and quaternionic two-qubit states. Then, we found those powers which fit our 
conjectured values (discussed in sec. II) for the associated Hilbert-Schmidt two-qubit sep- 
arability probabilities. The powers we found were 1.36743 {(3 = 1), 2.36904 (/3 = 2) and 
4.0632 ((3 = 4). We observe here a rough approximation to Dyson-index behavior, but 
can speculate that when and if the true forms of the Euler-angle separability functions are 
found, such behavior will be strictly adhered to. (In fact, one research strategy might be 
to seek functions that fully conform to this principle, while fitting the conjectured HS sep- 
arability probabilities. Also, below we will find closer adherence when we switch from the 
use of the participation ratio to a simple linear transform of the Verstraete-Audenaert-De 
Moor- function [69] (cf. [70]), which provides an improved bound on separability.) 

When we similarly sought to fit our prediction of ^ (that is, one-half of J| by the results 
of [62]) for the HS separability probability of generic minimally degenerate complex two- 
qubit states to a bivariate function proportional to a power of the participation ratio, we 
obtained a power of 6.11646, according rather well with Fig. 20. 

In terms of the Hilbert-Schmidt metric, the lower bound on the complex two-qubit sep- 
arability probability provided by the separable ball (-R(p) > 3) is rather small, that is 
23328V3 ~ 0-00272132, while relying upon the improved inequality reported in [69], 

VAD(p) = Ai - A 3 - < 0, (Ai > A 2 > A 3 > A 4 ), (68) 

it is 0.00365406. These figures are both considerably smaller than the comparable ones 
(0.3023 and 0.3270) given in [69] using (apparently) the measure (uniform on the simplex 
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of eigenvalues) first employed in [1]. 



E. Verstraete-Audenaert-De Moor function 



If we switch from the participation ratio R(p) to a simple linear transformation of the 
Verstraete-Audenaert-De Moor function, that is, 1 — VAD(p) (which varies over [0,1] for 
states outside the separable VAD set), in seeking to fit the trivariate separability function to 
our conjectured Hilbert-Schmidt two-qubit separability probabilities (sec. II), we find that in 

/ \ 3.15448 

the generic complex (f3 = 2) case, I 1 — VAD(p) J gives the best fit (for VAD(p) > 0). 

/ \ 1.53785 

Then, closely consistent with Dyson-index behavior, we obtained (1 — VAD(p)) as 
the best fit in the generic real {(3 = 1) scenario. (The VAD-bound (68) provides us with no 
useful information if we set A 4 = 0, so no relevance to the minimally degenerate two-qubit 
scenario is apparent.) 



F. Beta function fits to Euler-angle separability functions 

We can fit within 0.4% our conjectured values of ^ and ^ for the Hilbert-Schmidt 
separability probabilities of the complex and real minimally degenerate two-qubit states, 
respectively, by assuming-in line with the Dyson-index ansatz- that the Euler-angle separa- 
bility function in the real case is a regularized beta function (incomplete beta function ratio) 
[32, p. 11] of the form is( p )2(58, 22) (Fig. 21), and the Euler-angle separability function 
(Fig. 15) in the complex case, the square of that function. (Continuing along such lines, 
if we employ the fourth-power of the function, our estimate of the associated quaternionic 
separability probability is some 91.45% of the conjectured value of ^2^725 ~ 0.0386882.) 

Similarly, for the generic nondegenerate complex and real two-qubit states, we can 
achieve fits within 0.7% to both the conjectured HS separability probabilities of J| and 
Pp respectively, by taking in the real case the Euler-angle separability function to be 
If ^2(24,28) (Fig. 22) and its square in the complex case. (Use of its fourth power 

(l-VAD(p) 



to estimate the HS quaternionic two-qubit separability probability yielded a result 0.795969 

72442944 
936239725 



as large as the value, q 7 o 2 6 4 2 4 o 2 q 7 4 2 4 5 ~ 0.0795969, conjectured above (15).) 



40 



I. (58,22) 
^- t 

/ 0.8 

/ 0.6 

/ 0.4 

/ 0.2 

| 2 

0.7 0.75 0^8 0.85 0.9 0.95 ' S(p) 

FIG. 21: Incomplete beta function estimate of the bivariate Euler-angle separability function that 
closely reproduces the conjectured Hilbert-Schmidt minimally degenerate two-qubit real, complex 
and quaternionic separability probabilities 
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FIG. 22: Incomplete beta function estimate of the trivariate Euler-angle separability function that 
closely reproduces the conjectured Hilbert-Schmidt nondegenerate two-qubit real, complex and 
quaternionic separability probabilities 

IV. SUMMARY 

We have extended the findings and analyses of our two recent studies [11] and [17] by, 
first, obtaining numerical estimates of the separability function based on the (Euclidean, flat) 
Hilbert-Schmidt (HS) metric for the 27-dimensional convex set of quaternionic two-qubit 
systems (sec. II A). The estimated function closely conformed to our previously-formulated 
Dyson-index ((3 = 1,2,4) ansatz, dictating that the quaternionic {(3 = 4) separability func- 
tion should be exactly proportional to the square of the separability function for the 15- 
dimensional convex set of two-qubit complex ((3 = 2) systems, as well as the fourth power of 
the separability function for the 9-dimensional convex set of two-qubit real {j3 = 1) systems. 
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In particular, these additional analyses led us specifically to aver that 



sgUn) = (^) 2 ((3 - /i 2 )/i) 4 = (s™ pl M) 2 ,o < V < i 



(69) 




able to supplement (and fortify) our previous assertion that the HS separability probability 
of the two-qubit complex states is J| « 0.242424, claiming that its quaternionic counterpart 



cases (sec. II B), which still remain somewhat problematical. 

Further, we found (sec. II C) strong evidence of adherence to the Dyson-index ansatz 
for the 25-dimensional real and 35-dimensional complex qubit- qutrit systems with real HS 
separability function being proportional to a function of the form, 1 — (1 — v\v-i)i , with ui, u 2 
defined in (4). Subject to the validity of this separability function, we have obtained the cor- 
responding i? 2 constants ((5 — 1, . . . , 4) and estimated the complementary R 1 constants (the 
products -R1-R2 giving throughout, in our fundamental paradigm, the desired separability 
probabilities). 

Then (sec. II D), we determined that in terms of the Bures (minimal monotone) metric- 
for certain, basic simple scenarios (in which the diagonal entries of p are unrestricted, and 
one or two off-diagonal [real, complex or quaternionic] pairs of entries are nonzero)-that the 
Dyson-index power relations no longer strictly hold, but come remarkably close to doing so. 

Finally (sec. Ill), we examined the possibility of defining "separability functions" using 
the Euler-angle parameterization of Tilma, Byrd and Sudarshan [31], rather than the Bloore 
(correlation/off-diagonal scaling) framework [19]. Although now we are, prima facie, faced 
with a trivariate (in the complex two-qubit case) function, rather than a univariate one, 
we do not encounter the problem of having to determine an overall normalization factor 
for the separability function, since it is known that any density matrix with all eigenvalues 
equal to one another must be separable. It also appears that this simplifying feature further 
extends to the case where minimally degenerate (boundary) complex qubit-qubit states are 
considered (sec. Ill C), with any such state having its three non-zero eigenvalues all equal to 
I (lying on the boundary of the separable ball R(p) = 3) being necessarily separable. Use 
of the estimated Euler-angle trivariate separability function lent still further (numerical) 
support to the Jj conjecture [17] for the HS separability probability (and the associated J, 
conjecture for the hyperarea of the minimally degenerate two-qubit states) and the 168 °(^~ 1 ) 



is 



72442944 
936239725 



0.0773765. We have also commented on and analyzed the odd (3 = 1 and (3 = 3 
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"silver mean" conjecture [8] for the Bures separability probability of generic complex two- 
qubit states. 
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